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Abstract: We present numerical evidence for the existence of several types of static black 
hole solutions with a nonspherical event horizon topology in d > 6 spacetime dimensions. 
These asymptotically flat configurations are found for a specific metric ansatz and can be 
viewed as higher dimensional counterparts of the d = 5 static black rings, dirings and black 
Saturn. Similar to that case, they are supported against collapse by conical singularities. 
The issue of rotating generalizations of these solutions is also considered. 
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1. Introduction 

In recent years the interest in the properties of gravity in more than d = 4 dimensions has 
increased significantly. This interest was enhanced by the development of string theory, 
which requires a ten-dimensional spacetime, to be consistent from a quantum point of 
view. An unexpected result in this area was Emparan and Reall's discovery of the black 
ring in d = 5 spacetime dimensions [||, ^ . This asymptotically flat solution of the Einstein 
equations has a horizon with topology 5^ x S^, while the Myers-Perry black hole |^ has 
a horizon topology S^. This made clear that a number of well known results in d = 4 
gravity do not have a simple extension to higher dimensions. For example, the d = 5 
gravity allows for multi-black hole configurations regular outside and on the horizon. In 
this case, at least one of the constituents possesses a nonspherical topology of the horizon, 
the simplest examples being the black Saturn [Q] (a black ring with a central black hole), 
a diring ^ (two concentric coplanar black rings) and bicyling black rings Q (two black 
rings in orthogonal planes). 

However, while one can construct an encyclopedia of general relativity exact solutions 
in four and five dimensions, the situation for d > 5 is more patchy (see e.g. Q). For most of 
the cases, the known solutions are very special, with a large amount of symmetry. Moreover, 
it becomes clear that as the dimension increases, the phase structure of the solutions 
becomes increasingly intricate and diverse. The main obstacle stopping the progress in 
this field seems to be the absence of closed form solutions (apart from the Myers-Perry 
black holes), which were very useful in c? = 5. No general framework seems to exist for 
d > 5, and the issue of constructing black objects with a nonspherical horizon topology was 
considered by using various approximations or numerical methods. Most of the results in 
this area have been found by using the method of matched asymptotic expansions 10|. 



The central assumption is that some black objects, in certain ultra-spinning regimes, can 
be approximated by very thin black strings or branes curved into a given shape. However, 
this method has limitations; black holes whith no black membrane behavior (e.g. at high 



spins) would not be captured by this approach 

Although it would clearly be preferable to have analytic solutions^, some of the d > 5 
black holes with a nonspherical horizon topology can be constructed numerically, within 
a nonperturbative approach, as solutions of partial differential equations with suitable 
boundary conditions. 



^However, one should not exclude the possibility that most of these solutions will remain analytically 
intractable within a nonperturbative approach. 
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The main purpose of this paper is to present a general framework for a special class of 
static configurations with a symmetry group Rt x U{1) x S0{d—3) and to present numerical 
evidence for the existence of such solutions with nonspherical horizon topology. For d = 5, 
this framework reduces to that used in Q to construct generalized Weyl solutions. However, 
for higher values of the spacetime dimension, the solutions can be found only numerically. 
We argue that the basic properties of the d = 5 case still hold for d > 5 configurations with 
a symmetry group Rt x C/(l) x SO{d — 3), in particular the rod structure of the solutions. 
The simplest example of a d > 5 black object with a nonspherical horizon obtained within 



this approach was studied in Ref. |12| and has a horizon topology 5^ x S'^"'^. In this 
work, on the one hand, we extend these results and discuss the basic features of two new 
types of configurations representing composite black objects with (5"^ x 5°^"^) x S"^~^ - 
a generalized black Saturn - and also [S^ x S'^~^) x (5^ x 5"^~^) horizons -a generalized 
diring. On the other hand, within a slightly more general metric ansatz, we consider 
rotating solutions in either the 5^ or the S'^~^ spheres. 

This paper is organized as follows: in the next Section we present a systematic dis- 
cussion of this approach together with its limitations, while in Section III we present our 
numerical results. All solutions are found within a nonperturbative approach, by directly 
solving Einstein equations which for our ansatz reduce to a set of four nonlinear partial 
differential equations. 

Since all these solutions are plagued by conical singularities which seem to be unavoid- 
able in the absence of rotation, the issue of spinning solutions is addressed in Section IV. 
The results reported there are only partial, and so far we could not construct spinning 
balanced solutions. However, we expect that they will be useful for further work in this 
direction. 

Also, we have found that all new static solutions in this work have similar qualitative 
properties as their five dimensional counterparts. Therefore in the Appendix A we present 
the basic properties of the corresponding d = 5 solutions, which are known in closed form. 
Appendix B introduces a new coordinate system which simplifies the numerical calculations 
and leads to high accuracy (e.g. it has allowed to recover numerically the spinning balanced 
black ring starting with the static solution). 

2. The general formalism 

2.1 The field equations and a metric ansatz 

We consider the Einstein action 

/= — ^[ d'^x./^R d'^-^xV^K, (2.1) 



in a d— dimensional spacetime, with d > 5. The last term in ( p. ID is the Hawking-Gibbons 
surface term |T^, which is required in order to have a well-defined variational principle. K 
is the trace of the extrinsic curvature for the boundary dM and h is the induced metric of 
the boundary. Also, Gd is Newton's constant in d— dimensions; for simplicity, we shall set 
Gd = 1 in this work. 
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The upshot for our computations is that the hne element of the static solutions of 
interest can be cast in the following form (where < tp < 2it, < 6 < it/2 and dQ'^_^ the 
unit metric on S"^~^, while < p < oo, — oo < z < oo) 

The solutions constructed within this ansatz are static and axisymmetric, with a symmetry 
group Rt X U{1) X SO{d — 3) (where Rt denotes the time translation). While in principle 
it is possible to choose any kind of boundary conditions, we will only concentrate on black 
hole solutions which asymptote to flat spacetime. Moreover, the coordinates in (|2.2| ) have 
a rectangular boundary and thus are suitable for numerical calculations. 

A suitable combination of the Einstein equations, Gl = 0, Gp + = 0, = and 
= (with ip an angle on Qd-i), yields the following set of equations for the functions 
Uo, Ui, U2 

V^Uo + {VUo f + (Vf/o) • (VC/i) + {d- 4)(V[/o) • (VC/2) = 0, 

V^Ui + (VUif + (VC/o) • (VC/i) + (d - 4)(V[/i) • (VC/2) = 0, (2.3) 

V^U2 + {d- A){VU2? + (V[/o) • (VC/2) + {VUi) ■ (VC/2) -{d- 5)e2^-2^2 = q, 

and 

V^u - (VC/o) • (VC/i) -{d- 4){VUo) • (VC/2) - (fi - 4)(VC/i) • (VC/i) (2.4) 

+ i(d - 4)(d - 5) (e2--2^^ - (VC/2)2) = 0, 

for the metric function 1^, where we define 

(VC/) • (VV) = dpUdpV + d^Ud^V, V^C/ = d^U + 9^ C/. (2.5) 

The remaining Einstein equations Gz = 0, Gp — = yield two constraints. Following 
lH, we note that setting G^ = G^ = G^ + G^ = in V pG^'P = and V pG^''' = 0, we 
obtain the Cauchy-Riemann relations 

2dz (^/^G?) + dp {V^iGP - GD) = 0, 2dp {V^G^) - (V^(G^ - GD) = 0.(2.6) 

Thus, the weighted constraints satisfy Laplace equations, and the constraints are fulfilled 
when one of them is satisfied on the boundary and the other at a single point |14]. As we 
shall see, this is the case for all configurations discussed in the next Section. 

Although the Einstein equations take a simple form in terms of (Ui, z/), for the purposes 
of this paper it is more convenient to work with a set a functions fi defined as follows^ 

e^''^P''^ = h{p,z), e2^^(^'^) = /2(p,z), e^^^^P^^) = f,{p,z), e^^^^"'^) = /o(p, z). (2.7) 

This leads to a line element 

ds^ = -Mp, z)dt^ + fi{r, z){dp^ + dz^) + ^(p, z)dV'' + Hp, z)dnl_^, (2.8) 



^Some divergencies are avoided in this way. For example, fi — >■ would correspond to Ui — >■ —00 which 
is clearly not suitable for a numerical approach. 
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which was used in our numerical study of the d > 5 solutions. 

One might be concerned that ( |2.8| ) is too restrictive to leave room for new interesting 
black hole solutions. In higher dimensions, a priori this is not the case^, and, in the 
next Section, we shall present numerical evidence for the existence of nontrivial solutions 
which share the basic properties of some d = 5 configurations with a nonspherical horizon 
topology. 

Other more general metric proposals which may describe higher dimensional black hole 



solutions with a nonspherical horizon topology have been presented in [17|. However, due 
to their complexity, they will be out of our present scope. 

2.2 Known solutions 

2.2.1 Minkowski space-time 

In d > 5 dimensions, the flat spacetime metric can be written in the form (where < r < oo, 
0<V'<27r, 0<6l< 7r/2) 

ds'^ = -df + dr^ + r^idO'^ + sin^ Odip'^ + cos^ 9dnl_^), (2.9) 

thus with the metric on S'^~'^ written in terms of a warped product of S"^ and S"^~^. Then, 
for all dimensions, the coordinate transformation 

r = ^/2(p2 + ^2^1/4^ arctan(p/z), (2.10) 



leads to the equivalent form of (|2.9D 



1 



ds^ = -dt^ + ^ ^ ^ {dp^ + dz^) + ( + z) diP^ + ( Vp^+^ - z) dnl_^, (2.11) 
2^p^ + z^ 

where < p < oo, — oo < z < oo. 

An interesting observation here concerns the value of the determinant A for the non- 
conformal part of the metric^ {i.e. the line element (2.8) without the {p, z)-part). One can 
see that even for the simplest case of a Minkowski space-time within the parametrization 
(|2^) , A = -p^ for fi = 4, 5 only. 

2.2.2 Schwarschild-Tangherlini black hole 

The simplest example of a d > 5 nontrivial solution that can be studied within this approach 
corresponds to the Schwarzschild-Tangherlini black hole. Usually this metric is written in 
the form 

ds^ = -/(r) df + /(r)-^ dr^ + r^{de^ + sin^ 9 d^)'^ + cos^ 9 dnl_^), (2.12) 



^Black holes have to be of positive Yamabe type [|l^ and, if stationary, they have to be axisymmetric 

*The choice of this determinant has been proven to be crucial in recent progress on finding new classes 
of solutions ^] and also for the metric proposals in jl^] . 
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with /(r) = 1 — /i/r'^^^. This d— dimensional static black hole solution has an isometry 
group Rt X SO{d — 1). By a change of coordinates one can bring the metric to the desired 
conformal form ( |2.8D . The change of coordinates is 



which yields 



OL Oi 

p = -sin26' sinhG(r), z = - cos 26^ cosh G(r), (2.13) 



dp^ + dz'^ = y(cosh2G(r) - cos 40) i -^dr"^ + dO^ ] . (2.14) 



By simply integrating G(r) one finds 

G{r) = 2 [ ,[^dr = 2 [ , ^ dr = log[2^ (r^^ + y^r^-^ - ^i)^] + k, (2.15) 
J \ gee J ^Jf{r) r2 V J 

where k = —^133 log(4/i) is the integration constant. And, finally we fix a = 2^-3 fid-s to 
match asymptotically fiat space. 

The transformation ( p. 13 ) simplifies drastically in d = 5 

p=^-sin2e [l-^y^\\ z = i cos 20 (1 - ^) , (2.16) 
matching the findings of and in d = 7 where 

p = ^sm2e (1- ^y^^ r'^, z = ^cos2(9r^ (2.17) 

A straightforward but cumbersome computation leads to the following expression for the 
Schwarzschild-Tangherlini black hole in the (p, z) coordinates 

/ (d-3)/2 ^ Ud-3)/2^l)^ 

+ ^ 2.(.4l) + + '""^ + ^ 2.(.4l) + - 

where 

1/2 



y=l.(^p2 + Z^+P + V2^{p^ + z2)2 + ;x2(p2 _ ^2) + (p2 + ^2)^^ 

and P = V/"^ + - -z^) + (P^ + -z^)^- 

Moreover, one can also show that for d = 5 these expressions reduce to those in 

2.3 The rod structure of black hole solutions 
2.3.1 Five dimensional structure 



For d = 5, the coordinates in (|2.9D are the usual Weyl coordinates, while the sphere 
reduces to a single angular coordinate if, with < < 2tt. 
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In this case, it is most convenient to choose the three functions Ui as to satisfy the 
condition 



Y,U^ = log p. (2.19) 



This is compatible with the vacuum Einstein equations ( ^.31) , which for the choice ( 2.ig| ) 
imply also 

dp'^ p dp dz^ 

One can see that ( p.2Cl| ) is just Laplace's equation in a (fictions) three-dimensional fiat 
space with metric ds^ = dp^ + p^dQ^ + dz"^, whose solutions are well-known. 

From the other components of the Einstein equations Gp — = and Gp = 0, 
we obtain the equations which determine the function i'{p, z) for a given solution of the 
equation ( 2.20| ) 



i^' = -^^ + ^-{u[' + U',' + U'i-Ul-U^'-Ul), i^ = p{U[ + U', + U',\ (2.21) 

where a prime denotes the derivative with respect to p and a dot denotes the derivative 
with respect to z. Solutions with the ansatz (|2.2| ) and with C/i, U21 U3 and u satisfying the 
equations ( 2.20 ), (p.21|) are usually called generalized Weyl solutions p. 



This approach has proven very fruitful, a variety of physically interesting configura- 
tions being discussed in the literature. They can be uniquely characterized by the boundary 
conditions on the z— axis, known as the rod-structure ||2|, |l^, |l9|. One finds that the physi- 
cally relevant solutions for Ui can also be thought of as Newtonian potentials produced by 
thin rods of zero thickness with linear mass density 1/2, placed on the axis of symmetry 
in the auxiliary three-dimensional fiat space. Then the constraint (2.1£) states that these 
sources must add up to give an infinite rod. 

In this approach, the z— axis is divided into N intervals (called rods of the solution), 
[—00, zi], [zi,Z2],. . . , [zn-1, 00]. As proven in ||l^, in order to avoid curvature naked singu- 
larities at p = 0, it is a necessary condition that only one of the functions /o(0, z), /2(0, z), 
fsiOjz) becomes zero for a given rod, except for isolated points between the intervals. 

For the static case discussed here, a horizon corresponds to a timelike rod where 
/o(0,z) = while liuip^Q fQ{p, z) / p^ > 0. There are also spacelike rods corresponding to 
compact directions specified by the conditions fa{0,z) = 0, liuip^Q fa{p, z)/ p^ > 0, with 
a = 2, 3. A semi-infinite spacelike rod corresponds to an axis of rotation, the associated 
coordinate being a rotation angle. Demanding regularity of the solutions at p = imposes 
a periodicity 27r for both ■i/' and if. (However, when several ip- or 93-rods are present, it may 
be impossible to satisfy simultaneously all the periodicity conditions, see e.g. the examples 
in Appendix A). 

One of the main advantages of this approach is that the topology of the horizon is 
automatically imposed by the rod structure. This provides a simple way to construct a 
variety of solutions with nontrivial topology of the horizon (including multi-black objects). 
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bmce (|]2|) is linear, one can superpose different solutions for the same potential C/j. The 
nonlinear nature of the Einstein gravity manifests itself through the equation ( p.4| ) for the 
metric function v. 

2.3.2 Higher dimensional structure 

The central point in this approach^ is that the rod structure, as defined above for the d = 5 
case, can be used also for d > 5 solutions constructed within the ansatz ( |2.g| ). This fixes 
the boundary conditions along the z— axis for the functions fi and thus the topology of the 
horizon. 

However, note that the interpretation of a rod as corresponding to a zero thickness 
source with linear mass density 1/2, placed on the axis of symmetry in an auxiliary three- 
dimensional flat space is no longer valid for d > 5. Also, the relation ( ^.19 ) fails to be 



satisfied in this case, as one can see already for the simplest case of a Minkowski space- 
time. 

A crucial observation here is that, supposing the existence of a power series expansion 
in p, the Einstein equations imply^ the following form of the metric functions fi close to 
the z-axis, valid for any d > 5 

Mp, z) = fio{z) + p^fi2{z) + 0{p^), (2.22) 

where fikiz) are solutions of a complicated set of nonlinear second order ordinary differ- 
ential equations which we shall not present here. Then, similar to the d = 5 case, we 
suppose that the z— axis is divided into N intervals-the rods of the solution. Except for 
the isolated points between the rods, one assumes that only one of the functions fo{0,z), 
72(0, z), /3(0, z), becomes zero for a given rod, while the remaining functions stay finite at 
p = in general. (In fact, if more than one of these functions is going to zero for a given 
z inside a rod, one can prove following the arguments in [18| that we have a curvature 



singularity at that point.) Again, one imposes the condition that the intervals must add 
up to give an infinite rod. 

A finite timelike rod corresponds to an event horizon, where'' 

Up, z) = p2/o2(z) + p'Mz) + ..., h{p,z) = /io(z) + p^fuiz) + ..., (2.23) 

f2{p, Z) = f20{z) + p'f22{z) + . . . , h{p,z)= ho{z) + p' f32{z) + . . . • 

with limp_!.o p'^fi/fo = C3. As we shall see in the next Subsection, this fixes the Hawking 
temperature of the solutions. 

For a rod in the V'-direction, one finds the following expansion of the metric functions 
as /9 — )• 0: 

Mp, z) = fooiz) + p^fo2{z) + ..., Mp,z) = fioiz) + p'fuiz) + ..., (2.24) 

f2{p,z)=p'f22{z) + p'f24 {z) + . . . , hip, Z) = ho{z) + p^ h2{z) + .... 

''Some aspects of the proposal in this work can be found also in Ref. [^o[, which considers the numerical 
construction of the five-dimensional black rings with two independent angular momenta. 

"Such an expansion is also required by the regularity of the Kretschmann scalar at p = (but it does 
not guarantee it automatically). 



fik{z) here should not be confused with those in (2.24) 
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BH: S"*-2 horizon GBR: x S"^-^ horizon 

(a) (b) 



Q. ^^^^^ Q. 



-a a b 



GBS: (S'^ X S"*-*) X S"'-^ horizon GBD: {S^ x S^^-^) ^ (5-2 x S'^-'^) horizon 

(c) (d) 



£1 n 



Figure 1. Rod structure of several static black objects in d > 5 spacetime dimensions. These 
include the Schwarzschild-Tangherlini black hole (BH), the generalized black ring (GBR), the gen- 
eralized black Saturn (GBS) and the generalized black di-ring (GBD). In the diagrams the thin 
lines represent the z— axis and the thick lines denote the rods. 

The important feature here is that the constraint equation Gp = implies fio{z)/f22{z) = 
ci > 0, i.e. a well-defined periodicity for the coordinate ip. The value of ci is not fixed 
apriori and follows from the details of the solutions. 

For d = 5, a similar result is found when interchanging /2 and /a, i.e. for a rod in the 
(/^-direction. The periodicity of 99 there is arbitrary, being again fixed by the constraint 
equation = 0, i.e. limp_j.o /o^/i/Zs = C2 > 0. However, the picture is very different^ for 
d > 5, in which case we find 

h{p, z) = foo{z) + p^fo2{z) + ..., fi{p,z) = fioiz) + p^fuiz) + ..., (2.25) 

f2{p, Z) = f20{z) + p^f22{z) + p^f2A{z) + . . . , h{p, z) = p^ f^iz) + ..., 

i.e. the Einstein equations impose the following requirement for an il— rod: limp^g p'^fi/fs = 
1, which is an important new feature. As we shall argue in the Section 2.5, this condition 



Formally, this is a consequence of the presence of the (d — 5) coefficients in the eqs. (2.3), (2.4) 
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prevents us to construct d > 5 static black rings (or multi-Schwarzschild-Tangherlini black 
holes) within the ansatz ( |2.8D . 

Thus, depending on the physical situation we consider, the boundary conditions along 
the z— axis are fixed by the above relations. The obvious boundary conditions for large 
p, \ z\ are that fi approach the Minkowski background functions (which are read from 2.11j ): 



1 



/o(p,z) = l, /i(p,z) = ^==, f2ip,z) = ^/p^ + z^ + z, /3(p,z) = vV+^-^-(2.26) 

This is in fact the simplest solution of the Einstein equations in d > 5 dimensions, as we 
have seen already in the Section 2.2.1. There the function /2 vanishes for p = 0,z < 0, 
and /s for p = 0, z > 0, which, in the language of the Weyl formalism, corresponds to two 
semi-infinite rods [— oo,0] and [0, (X)]. 

Similar to the d = 5 case, the topology of an event horizon is fixed by the boundary 
conditions satisfied by /2 and at the ends of the corresponding (finite) timelike rod. 
For example, if either end of this rod continues with rods of different directions {^^J and 
O), then the event horizon has an S'^~'^ topology and (for a single horizon) the solution 
is a Schwarzschild-Tangherlini black hole (see Figure la). A black object with 5^ x S"^~^ 
topology of the horizon is a 'generalized black ring' and has the metric function /2 vanishing 
at both ends of the finite timelike rod associated with the horizon, see Figure lb. (For d = 5, 
this corresponds to the static black ring in (P.) One can consider as well a 'generalized 
black Saturn' combining both types of black objects above, with two different horizons 
(Figure Ic). The rod structure for a solution consisting of two black holes, both of them 
with S'^ X S"^~^ topology of the horizon (thus a 'generalized black diring') is shown in Figure 
Id. The basic properties of these objects are discussed in the next Section. 



It is tempting to conjecture that, similar to the d = 5 case |19], a d > 5 solution within 
the ansatz (2.8), is uniquely specified by its rod structure. 



2.3.3 Physical quantities 

The ADM mass M of the solutions^ can be read from the asymptotic expression for the 
metric function gu 

levrAf 

id- 2)y,_2(p2 + .2)(.-3)/4 + • • • • (2-27) 

Supposing we have an event horizon at p = for some zi < z < Z2, the horizon metric is 
given by^'' 

da^ = MO,z)dz^ + f2{0,z)dil^^ + f3{0,z)dnl_^. (2.28) 

Two quantities associated with the event horizon are the event horizon area Ah and the 
Hawking temperature. For the metric ansatz (|2.8| ) these are given by 



Ah = AVVd_4 dzJh{0,z)h{0,z)ftH0,z), Th = ^ lim J {°/^'^\ , (2.29) 
Jzi ^ 2tt P^oV p^fi{p,z) 



®The discussion here follows the general formalism in which contains also several applications for 
d = 4, 5 exact solutions. 

^"if there are several horizons, then one should write such an expansion for each of them. 
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where Vd-4 is the area of the unit sphere S'^~'^ and A-^ the periodicity of the angular 
coordinate '0 on the horizon. 

A solution with n different event horizons satisfies the Smarr law 

^=4|!rl)X^^i''^2^ (2.30) 

Considering now the case of a space-like ■;/;— rod for some < z < Z4^, one starts by 
writing the line element on this {d — 2)-dimensional surface S 

da^ = MO, z)dz^ + /3(0, z)df]l_4 - /o(0, z)dt'. (2.31) 

The first quantity of interest is the proper length of the rod 

Z4 



L = / dzV/i(0,z), (2.32) 



which, for a finite rod, differs from the coordinate distance = 2:4 — ^3 (although it is 
proportional to it). 

All solutions in this work possess a conical singularity for (at least) a region of the 
z-axis. To define a conical singularity for a rotational axis with angle if), one computes the 
proper circumference C around the axis and its proper radius R and defines 



dC 
''=dR 



= hm #^ = hm ^^^^^ (2.33) 



where l^ip is the period of ip. The asymptotic flatness implies l^ip = 2tt. Then the presence 
of a conical singularity is now expressed^^ by means of: 



^ = 27r - a = 27r ( 1 - lim J {'/^' ^\ ] , (2.34) 

such that 6 > corresponds to a conical deficit, while S < corresponds to a conical excess. 
A conical deficit can be interpreted as a 'string' stretched along on a certain segment of 
the z— axis, while a conical excess is a 'strut' pushing apart the rods connected to that 
segment (in fact, the 'struts' and 'strings' are (d — 3)-dimensional surfaces, i.e. higher 
dimensional analogues of the d = 4 cosmic strings). Also, a constant rescaling of ip can 
be used to eliminate possible conical singularities on a given segment, but in general, once 
this is fixed, there will remain conical singularities at other -(/^-segments. For all solutions 
in this work, we have prefered to set the conical singularity on a finite ip-rod such that our 
solutions are asymptotically flat, meaning that Aip = 2tt. Moreover, in the presence of a 
conical singularity, the manifold M naturally factorizes as = Cq x S, where Ca is the 
two-dimensional conical surface p — tp and S is the remaining {d — 2)-dimensional surface, 
which may be seen as the world- volume of the conical defect p2|. 



^^Note that, in some sence, fixing S is tlie analogue of computing the Hawking temperature on the 
Euchdean section. 
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For practical reasons, we have found it convenient to introduce the quantity 

e_ V(27r) 
1 - V(2vr) ' 



(2.35) 



which has a finite range (5 — t- — 1 for 5 — )■ ±00) and measures the 'relative angular 
deficit /excess'. 

As argued in 23], the asymptotically flat black objects with conical singularites still 



admit a thermodynamical description. Moreover, when working with the appropriate set 
of thermodynamical variables, the Bekenstein-Hawking law still holds for such solutions. 
The mass-energy which enters the first law of thermodynamics does not, however, coincide 
with the ADM mass; it differs from the latter by the energy associated with the conical 
singularity, as seen by an asymptotic, static observer. 

The Ricci scalar of M, ("^i?, can be represented near S in the following local form 

Hi? = i? + 2(2^-a)5s, (2.36) 

where R is the curvature computed in the standard way on the smooth domain Ai/Ti of 
^A. Here, (5e is the Dirac 5-function, with /^e = /s /• ^ direct integration of (|2.36| ) 



leads to [g4| 



[ (")i?= [ R + 2{2tt - a)Area, (2.37) 

Jm Jm/t. 

where Area is the area of S, i.e. the space-time area of the conical singularity's world- 
volume. For the metric ansatz in this work, the expression of the extensive parameter 
Area is 

/3Vd_4 r dz^fo{0,z)MO,z)ft\0,z), (2.38) 



Area 

(5 = 1/Th being the periodicity of the Euclidean time. 

In the presence of conical singularities, the solutions cannot be viewed as vacuum 
solutions and there is a matter source (the struts) which supports the conical singularities. 
The stress energy tensor associated with the singularity can be computed by using the 



Einstein equations Gij = SnTij. The results in |22| also show that the singular part of the 
Ricci tensor has components only in the p — ip plane, such that i?^ = for the remaining 
components. It follows that the only non-vanishing components of T/ are 

^' = -'^''^(2^-«Ms, with {i,j)^{p,y,). (2.39) 



A direct consequence of this result is that the conical deficit /excess as defined by ( 2.34| ) 



6 = 2Tr — a, corresponds to the pressure exerted by the strut. This is found by integrating 
the TJ^-component over Ca 



P 
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Moreover, as seen from (|]3|), the energy density /i of the matter source supporting the 
conical singularity is also // = —6/8tt. Thus 6 < (the case of the solutions in this work) 
corresponds to a negative energy density source. 

Another quantity of interest is the total energy associated with the strut as seen by 
a static observer placed at infinity. This is found by integrating the T/-component over a 
t = Iq = constant hyper-surface, 



where we have defined 



Area , , , 

A^—^. (2.42) 



Thermodynamics of a system with a conical singularity in the bulk can be approached 
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by using the path-integral formulation of quantum gravity [13|. The first step here is to 
evaluate the total tree level Euclidean action of the system. The new feature introduced 
by the conical singularity is to add an extra contribution to this quantity which can be 
evaluated by using the relation ( p.37| ). Then the total action is 

^ = ^0 - |^^/3, (2.43) 



where Iq is the usual tree level action [13| found when neglecting the conical singularity. 

As argued in [^], the first law of thermodynamics for static vacuum solutions with a 
conical singularity reads 

dM = TudS + TdA, (2.44) 



where A (as defined by (2.42)) is the extensive parameter which takes into account the 
presence of the conical singularity and T = P = —6/8tt is the associated "tension". 

In a canonical ensemble, one keeps the Hawking temperature Th and the extensive 
parameter A fixed. The free energy F is 

F[Th, A]=ThI = M- ThS. (2.45) 

Then the entropy 5, mass M and tension T of the physical system are given by 

(2.46) 



S 



OF 

, M = F + ThS, T= — 
A 9A 



Th 



This approach has been applied in for several d = 4, 5 static solutions with conical 
singularities in the bulk which are known in closed form, the effects of rotation being 
considered in [^]. The results there show that S = Ah/^ in all cases, as expected. 

In principle, one can use the same approach to discuss the thermodynamics of the 
solutions in this work. The only complication appears for multi-black objects. For thermal 
equilibrium, the individual black holes should have the same Hawking temperature. More- 
over, if there are several different finite V'-rods, then the corresponding tension parameters 
5i should be equal. 

^^In fact, it is possible to define tlie thermodynamic quantities also with the conical singularity stretching 
towards the boundary, see e.g. |25[. 
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2.4 Numerical procedure 

All new solutions in this work are found within a nonperturbative approach by solving 
numerically a set of four nonlinear partial differential equations. These equations have 
dependence on two variables and are subject to suitable boundary conditions which follow 
from the required rod structure and asymptotic flatness. 

In this scheme, the input parameters provided to the solver are the positions of the 
rods and the value d of the spacetime dimension. The numerical integration eventually 
converges and provides an output consisting of the functions fi and their first and second 
derivatives with respect to p and z. The relevant physical data {e.g. the ADM mass, the 
Hawking temperature(s) etc.) are computed from this numerical output. 

In practice, we have found it convenient to take 

f^ = f^Fi, (2.47) 

where are some background functions^^, given by the metric functions of a d = 5 
solution with some rod structure. The numerical computation is performed working with 
the functions Fi. The advantage of this approach is that, since Fi > everywhere, the 
functions fi will automatically satisfy the desired rod structure also for the d > 5 solutions. 
Moreover, this choice 'absorbes' the divergencies of the functions /2 and as p ^ oo, 
z —7- iboo originating in the imposed asymptotic behaviour (and also the divergencies of /i 
at the end of the finite rods^^). 

The equations satisfied by Fi can easily be derived from those satisfied by As for 
the boundary conditions, the relations ( ^.24 ), ( 2.23| ) imply 



dpFi\p=Q — 0, 

for a V'-rod or an event horizon, and 

apFo|p=o = apFi|p=o = 5pF2|p=o = 0, Fi|p=o = i^3|p=o, (2.48) 

for an il-rod. As /3 — )• oo or z — )• itcxD one imposes the obvious condition Fi = 1. 

The constraint equation Gp = implies F2/F1 = const, on a ip-iod. Now, to be 
consistent with the assumption of asymptotic flatness, one finds const. = 1 for a ^— rod 
extending to infinity. The value of this ratio for a finite rod results from the numerical 
results. 

In the numerical calculations, one starts with the corresponding d = 5 solution, (i.e. 
Fi = 1), and slowly increases the parameter d. This leads to nontrivial deformations of Fi. 
The physical values of d are integers. For all types of black objects in this work, we have 
studied solutions in d = 6 dimensions in a systematic way. A number of d = 7 solutions 
were also obtained but we did not fully investigate this case except for single black objects. 



Although sometimes it is not stated expUcitly, the use of a suitable set of background functions is a 
common feature of numerical studies in general relativity. For example, for axisymmetric configurations, 
this factorizes the trivial angular dependence of the metric functions. 

""^^Note that for all d = 5 solutions the function /i(0, z) behaves as l/\z — Zi\ as 2 — > Zi, with Zi the end 
point of a finite rod. This behaviour is recovered for all d > 5 solutions in this work. 
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We think that solutions with d > 7 are also very likely to exist; however, we could not 
construct them within the approach in this work and their study may require a different 
numerical method^^. 

In practice, one introduces also new compactified coordinates x = p/{l + p), u = 
arctan(z), with < x < 1, — 7r/2 < u < tt/2. The equations for Fi are then discretized 
on a non-equidistant grid in x and u. Typical grids used have sizes around 80 x 150, 



covering the full integration region. (See p9| | and [30| for further details and examples for 
the numerical procedure.) All numerical calculations are performed by using the program 
FIDISOL/CADSOL ||2^ , which uses a Newton-Raphson method^^. This software provides 
also an error estimate for each unknown function. For d = 6, the typical numerical error 
for the functions is estimated to be lower than 10~^. This error increases to several percent 
for most of the d = 7 solutions. 

Perhaps the most problematic aspect of this approach is that the derivatives of the 
functions Fi diverge at the end points of the rods (although the functions Fi are smooth 
everywhere). However, this appears already for the d = 5 exact solutions^^ and is an 
unavoidable feature of the Weyl-type coordinate system ( |2.8D . To assure that these diver- 
gencies are coordinate artifacts, we have verified that the Kretschmann scalar stays finite 
everywhere, in particular at p = (here we ignore the 5-Dirac terms in the expression of 
Riemann tensor in the presence of a conical singularity ||2^). However, this non-smooth 
behaviour makes the numerics more involved as compared to other problems which were 
solved with similar methods and the same software. In particular, this approach requires 
a careful construction of the mesh in the region close to the end points of the rods. One 
should also remark that the presence of a conical singularity for the solutions in this work 
has a rather neutral effect on the numerics, since the solver does not notice it directly. 

Furthermore, as a test of the software, we have recovered numerically the d = 6, 7, 8 
exact Schwarzschild-Tangherlini black hole starting with the d = 5 solution, which, for 
the ansatz ( ^^ ) possesses already the basic features of other more complex black objects. 
Another accuracy test of our solutions was provided by the Smarr relation ( p. 30] ) . 

A further numerical test is presented in Appendix B, where the d = 5 balanced black 
ring is recovered numerically. The background functions in this case correspond to those 
of the static exact solution. 



We believe this is not a worrisome aspect. For example, a similar situation was found in the past for 
nonuniform black strings in Kaluza-Klein theory, which at the beginning could be constructed numerically 
for d = 6 only Jl^. However, subsequent work has managed to extend these solutions to d = 5 and 
7 < d < 11 [^^. Kaluza-Klein caged black holes were so far constructed within a nonperturbative approach 
for d= 5,6 only [|||. 

^^Thus providing a good initial guess solution is essential for the convergence of the numerical process. 

^'^For a d = 5 black ring, the derivatives of the metric functions are finite everywhere when writing 
the solution in a ring-coordinate system. Although one can devise a ring-coordinate system (or more 
complicated versions adapted to the solutions in this work) for d > 5 as well, we could not use it in practice 
within a numerical scheme. The main problem seems to be that the spacelike infinity is approached for a 
single point in that case. 



-15- 



2.5 The issue of d > 5 black ring solutions 

One may expect that the approach proposed in this work can be used to construct a static 
black ring as well. But to be able to do so the metric has to have an unbroken SO{d — 3) 
symmetry group. The approximate solutions for thin black rings Q do indeed preserve 
this symmetry and are cohomogenity-2. Therefore, we expect their static metrics in the 
(/O, z) coordinates to be within our ansatz. 

The rod structure of these solutions in the generalized approach we are considering 
would consist of a semi-infinite space-like fi— rod [— oo, — a] (with /3(0,z) = there), a 
finite time-like rod [—a, a] (the horizon, /o(0,^;) = 0), a second (finite) J7-rod [a,h] and a 
semi-infinite space-like -i/^-rod [6, oo] (with a vanishing /2(0, z)). Thus, it is obvious that 
the horizon topology of such an object would be S'^~^ x . 

In principle, one may consider as well multi-black objects with several fi— rods, the 
simplest case corresponding to a two Schwarzschild-Tangherlini black hole system, both 
with a 5"^~^ topology of the horizon. Indeed, such solutions are known in closed form in 



four ||3l|, and five dimensions |33|] . 



On general grounds, we know that higher dimensional black rings are possible and 
Ref. 1^ found an approximate construction of such solutions based on the matched asymp- 
totic expansion method. We note that the solutions in Ref. describe balanced spinning 
black rings, but unbalanced and, in particular, static solutions should also exist. The 
analysis in Ref. Q shows, that when the equilibrium condition is not satisfied, naked sin- 
gularities arise in the plane of the ring. For higher-dimensional rings, the singularities 
are thus necessarily stronger than in five dimensions, where conical singularities suffice to 
balance the rings. 

This observation is in correspondence with our analysis in ( |2.25| ), which shows that for 
d > 5 one cannot assign to ^d-i ^'^J higher dimensional counterpart of a conical defect 
since limp_).o p^/s/Zi = 1 in that case. One may conclude that the static black rings do 
not admit at p = a power series expansion of the form ( 2. 221 ), since the static limit of a 



higher dimensional back ring will possess a stronger singularity than a conical one. 

Clearly, for the envisaged numerical construction of such higher-dimensional static 
black rings the presence of curvature singularities provides a strong obstacle. Therefore 
we conclude that we cannot study d— dimensional static black rings within the metric 
in a straightforward manner and a different approach seems to be required. 



because 



^*Moreover, a d— dimensional static black ring might not be described within the ansatz (2.1 
the sphere Q.d-i could be deformed in that case (i.e. the metric functions would have a dependence on (at 
least) one coordinate on that sphere). 
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1/f, 





Figure 2. The profiles of the metric functions for a d = 8 Schwarzschild-Tangherhni solution 
with a = 1. 



3. General results for static black holes 
3.1 Uni-horizon black holes 

3.1.1 Schwarzschild-Tangherhni black hole: S'^^^ horizon 

The simplest example of a d > 5 black object that can be studied within this approach 



corresponds to the Schwarzschild-Tangherhni black hole. In Section 2.2.2 we have shown 
its form within the metric ansatz (|2.8| ). 

The rod structure of this black hole consists of (see Figure 1) 

• A semi-infinite space-like rod [— oo,— a] in the d/dip direction (with /2(0,z) = 
there) , 

• A finite time-like rod [—a, a] in the d/dt direction (the horizon, /o(0,z) = 0), 

• A semi-infinite space-like rod [a, oo] in the fi-direction (with a vanishing /3(0,z)). 
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Thus the topology of the horizon is S'^~'^ as required. Requiring the absence of a conical 
singularity imposes a periodicity 2tt for the coordinate ip. The only parameter here is 
a > 0, which fixes the ADM mass of solutions (the only global charge for a Schwarzschild- 
Tangherlini black hole). 

To get the feeling for the type of functions we will be dealing with later in the numerical 
computations, we show in Figure 2 the form of the metric components, namely the metric 
functions /j, exhibiting this rod structure for a typical d = 8 Schwarzschild-Tangherlini 
solution. 

In principle, most of the physically relevant properties of the Schwarzschild-Tangherlini 



black hole can also be rederived within the metric ansatz (2^). However, the required 



computation is much more difficult for that coordinate system. 
3.1.2 Generalized black ring: S"^ x 5"^ horizon 

Very likely, these are the simplest d > 5 black objects with a nonspherical topology of the 
horizon which can be constructed within a nonperturbative approach. 
The rod structure in this case consists of (see Figure lb) 

• A semi-infinite space-like rod [— oo,— a] in the d/dip direction (with 72(0,2;) = 
there) , 

• A finite time-like rod [—a, a] in the d/dt direction (the horizon, /o(0,z) = 0), 

• A a second (finite) ip-rod [a,b] (with /2(0,z) = 0), 

• A semi-infinite space-like rod [6,00] in the il-direction (with a vanishing /3(0, z)). 

The problem has two length scales a and b, roughly corresponding to the event horizon 
radius and the radius of the round r2rf_4-sphere. These two input parameters fix the ADM 
mass M and the tension 5 which are the global charges here. 

The basic properties of these solutions were discussed in for d = 6,7. For d = 5, 
they correspond to the static black ring found in (its basic properties are reviewed in 
Appendix A.l). The d > 5 configurations share the basic properties of the five dimensional 
counterparts. This is why we propose to call them 'generalized black rings'. 

Unfortunately, all static solutions with a x S'^~^ horizon topology have a conical 
excess 6 on the finite ip—rod. In terms of the dimensionless ratio a/b, the generalized black 
rings smoothly interpolate between two limits. First, as a/b — )• 1, one finds 6 — )• —1 {i.e. the 
conical excess 5 — )• —00) and the Schwarzschild-Tangherlini metric is approached (the finite 
■i/'-rod vanishes). As the second ^z^— rod extends to infinity {a/b — )■ 0), the radius on the 
horizon of the round 5"^~^-sphere increases and asymptotically it becomes a. {d — 4)-plane, 
while J — )• 0. After a suitable rescaling, one finds the four dimensional Schwarzschild black 
hole uplifted to d dimensions {i.e. a black {d — 4)-brane). 
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Figure 4. Features of generalized black ring solutions, (a): The scale free ratio T^^A^^j 
is shown as a function of the relative angular excess d — S/{2tt — S). The parameter ad = (d — 
3)2(5-2d)/(d-2)^(3-d)/2(d-2)p(^^ _ '^y2y/{2-d) j^^^g ^gg^^ choscu such that the point (1, -1) on the 

plot corresponds to the Schwarzschild-Tangherlini black hole, (b): The scale free ratio A/A4 is 
shown as a function of the ratio between the two length scales b/a. 

The Kretschmann scalar Rij^^iRJ'^^^ of typical d = 5, 6 solutions is shown in Figure 3. 
There one can see the nonsingular character of this type of configurations (here we do not 
consider the (5-Dirac terms in the expression of Riemann tensor in the presence of a conical 
singularity (2^). 

In Figure 4 we show the scale free ratio Tjj^A^Jj^'^ as a function of the relative angular 
excess 5 = 5/{2-k — 5) for d = 5, 6, 7 solutions. The dimensionless ratio A/M. is also shown 
there as a function of the ratio h/a of the two length scales. One can see that, as expected, 
the pattern for d = 5 repeats for the higher dimensional configurations. Further details on 
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these solutions including typical profiles of the functions , Fi are given in Ref . 
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3.2 Mult i- horizon black holes 

The general ansatz proposed in Section 2 offers the possibility to construct multi-black hole 
solutions as well. As discussed above, a limitation of this approach is that all configurations 
would have a single fi-rod extending to infinity. (For example, this prevents us from 
constructing d > 5 multi-Schwarzschild-Tangherlini solutions.) However, the number of 
rods in the -0 or t-directions are not constrained (although the numerical accuracy decreases 
with the number of rods). 

In what follows, we present numerical evidence for the simplest asymptotically flat 
multi black objects, corresponding to 'generalized black Saturns' and 'generalized black 
dirings' (although more complex configurations are very likely to exist). 

3.2.1 Generalized black Saturn: (5^ x S'^~'^) x S'^'^ horizon 

The d = 5 static black Saturn describes a multi black hole configuration, with a black 
ring with horizon topology S"^ x around a Schwarzschild-Tangherlini black hole. This 
exact solution is found as a static limit of a rotating solution originally presented in Q (see 
Appendix A. 2 for a review of its basic properties). 

We are interested in the generalization of this type of configuration to d > 5 within 
the metric ansatz ( p.Sp . There we have again two black objects, with a topology of the 
horizon 5^ x S"^~^ for the generalized black ring and S"^~^ for the central black hole. (Note 
that the horizon of the central black hole is not a round sphere). 

This type of solution is found by imposing the following rod structure: 

• A semi-infinite space-like rod [— oo, — a] in the d/dip direction (with /2(0,z) = 0), 

• A finite time-like rod [—a, a] in the d/dt direction (/o(0,2;) = 0), 

• A finite space-like rod [a,b] in the d/dip direction (where /2(0,z) = again), 

• A second finite time-like rod [6, c] in the d/dt direction, (/o(0, z) = 0), 

• A semi-infinite space-like rod [c, oo] in the 0-direction (with a vanishing /3(0, z)). 

This rod structure is illustrated in Figure Ic. 

The profiles of a tyical d = 7 solution are shown in Figure 5. Thus the problem has 
three input parameters a, b, c fixing the mass of the individual components and the distance 
between them. In practice, one can always fix one of the parameters a, b, c {i.e. the length 
scale of the problem) and vary the other two. 



-20- 




- 21 - 



I 1 1 1 1 1 -0.1 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 2.2 2.4 




Figure 6. d = 6 generalized black Saturn solutions: a number of relevant dimensionless quantities 
are shown as a function of the ratio b/c for several fixed values of a/c (left) and for a varying ratio 
c/b for fixed a/b (right). In these plots the indices BR, BH stand for the objects with S"^ x 5''*^'* 
and S'^^'^ horizon topology, respectively. 

As expected, all our solutions have a conical singularity on the rod between the hori- 
zons, with (5 < there, which corresponds to a conical excess (and thus a negative energy 
density for the strut source). 
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Starting with the dependence of solutions on the ratio h/c for fixed a/c one can see that 
the solutions smoothly interpolate between two limits. First, for h/c ^ 1, the horizon with 
S'^~'^ topology vanishes and the solution reduces to the higher dimensional generalizations 
12] of the static black ring in P]. Another limit of interest is 6/c — t- a/c, in which case the 



finite ip-Tod vanishes and a Schwarzschild-Tangherlini configuration is recovered. 

We have studied as well the dependence of the solutions on the ratio c/b for fixed a/h. 
There, for c/b — )• 1, the horizon with S"^~^ topology vanishes and the solution reduces to 
a generalized black ring. Other interesting limits are a — ?• 5 or a = 0, in which cases the 
Schwarzschild-Tangherlini configuration is recovered. 

Some results of the numerical integration supporting the above statements are shown 
in Figure 6 for several fixed values of a/c and a varying h/c (left) and for a varying ratio c/h 
for fixed a/h (right). One can see that the generic solutions are not in thermal equilibrium, 
7^ "^H^ (also, we could not find configurations with a vanishing Hawking temperature 
for one of the components). For all solutions, we have noticed a good qualitative agreement 
of their behaviour with that found for d = 5 static black Saturns. 

3.2.2 Generalized di-rings: {S^ x S"^^) x {S^ x S'^-^) horizon 

For d = 5, the Einstein equations have an exact solution describing two concentric black 
rings In the static limit, this asymptotically flat configuration is supported by a 

strut with positive pressure and negative energy density. The explicit form of the d = 5 
di-ring solution is given in Appendix A. 3. 

By using the same approach as in the previous cases, we could construct higher dimen- 
sional generalizations of this static configuration. There we have again two black objects, 
both with a topology of the horizon S"^ x S'^~'^. They are found for the following rod 
structure: 

• A semi-infinite space-like rod [—oo, —d] in the d/dtp direction (with /2(0, z) = 0), 

• A finite time-like rod [— d, — c] in the d/dt direction (/o(0, z) = 0), 

• A finite time-like rod [— c, —a] in the d/d^ijj direction (where /2(0, z) = again), 

• A second finite time-like rod [—a, a] in the d/dt direction (/o(0,z) = 0), 

• A finite space-like rod [a, 6] in the d/d^ij: direction (/2(0, z) = 0), 

• A semi- infinite space-like rod [6, oo] in the il-direction (with a vanishing /3(0, z)). 

This rod structure is illustrated in Figure Id. 

Given the presence of four finite rods, finding such solutions is a more difficult problem, 
and we did not manage to obtain d = 7 numerical solutions with reasonable accuracy. 
However, we think this is due to the limitations imposed by our approach only. 

The metric functions fi and the functions Fi used in the numerical calculations change 
smoothly with the rod parameters a, b, c and d. Typical profiles of the solutions are 
presented in Figure 7 as a function of z for several values of p. 
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Figure 8.^ = 6 generalized black; diring solutions: a number of relevant quantities are shown as 
a function of the ratio a/b for several fixed values of c/6, d/c (left) and for a varying ratio b/a for 
fixed a/c, d/c (right). In these plots the indices 1, 2 stand for the objects with an event horizon at 
p = 0, — d < z < — c and p ~ 0, — a < z < a, respectively. 

One can see that the functions Fi are smooth outside of the z— axis and show no sign 
of a singular behaviour (although they have a complicated behaviour at p = 0). 

In this case we have studied the dependence of the solutions on the parameters a, b, c, d 
for two different situations (In fact a simple rescaling leads to a dependence on only three 
dimensionless quantities.). In the first case (see Figure 8 (left)), we have fixed b,c,d and 
varied the parameter a associated with the length of the second timelike rod (one of the 
constants b, c, d can be taken to represent the length scale of the problem). 

As a —7- 0, a generalized black ring is approached since the second horizon vanishes, 

(2) (2) 

which implies — )■ oo and Aj^ — t- 0. As a/c —t- 1 for fixed b,c,d, the first finite Tp-rod 
vanishes and the two horizons coalesce to form a single black object with S"^ x S'^~^ horizon 
topology. 

The second case we have investigated corresponds to fixed a, c, d and a varying b (see 
Figure 8 (right)). The relevant limits here are 6/a — >• 1 and b/a ^ oo. As a/6 —t- 1, one finds 
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62 — —1 (i.e. the conical excess 62 — )■ —00) and a generalized black Saturn configuration 
is approached (the second finite ^— rod vanishes). The case b/a — )• 00 is more involved 
and we have found it more difficult to investigate this limit. As the second ■0— rod extends 
to infinity, one expects to recover, after a suitable rescaling, the four dimensional double- 
Schwarzschild configuration (i.e. the Bach-Weyl solution |3^]) uplifted to six dimensions. 
This configuration still has a conical singularity in between the two black holes which 
provides the force balance that allows its existence. 

Again, we have noticed a good qualitative agreement of this behaviour with that found 
for the d = 5 static dirings exact solution (see Appendix A3). The generic solutions have 
different Hawking temperatures for the individual components, T^^ ^ '^h^ ^ ^'^'^ thus are 
not in thermal equilibrium. Also, it seems that there are no generalized black diring 
solutions with a vanishing Hawking temperature for one of the components. 

4. Remarks on rotating black holes with a nonspherical horizon topology 

All static solutions with a nonspherical horizon topology discussed above are plagued by 
conical singularities, which seems to be an unavoidable feature of all such asymptotically 
flat black objects. 

However, for their d = 5 counterparts, the conical singularities are eliminated by 
spinning the configurations, the rotation providing the centrifugal repulsion that allows a 
regular solution to exist. 

On general grounds, one expects the d > 5 new solutions in this work to possess 
rotating generalizations. Thus one may hope that by adding (at least) one spin to the 
system (without changing the rod structure) the configuration will be balanced {i.e. without 
conical singularities) for a critical value of the angular momentum. Unfortunately, there is 
no simple procedure to spin a given static system. Moreover, all known techniques fail for 
the solutions considered in this work. Therefore, we will again use a numerical approach. 

In what follows, we present two different proposals for a metric ansatz which may 
describe rotating black holes with a nonspherical horizon topology. The first proposal 



is a straightforward extension of the ansatz (2.8) used in the static case, and leads to 
equations with dependence on two variables. The second ansatz proposal implies an a 
priori dependence of the unknown metric functions on three variables. 

4.1 A generalized metric ansatz 

As discussed in [^], the Weyl coordinates and the rod-structure employed to construct 
d = 5 static axisymmetric solutions can be generalized to the rotating case. Here we 
present arguments that the general framework proposed in Section 2 can also be generalized 
to include rotation. 

In the simplest case, one considers a slightly more general metric form than (|2.2| ), with 

ds^ = e^^^f^'^dp^ + dz^) + Gijip, z) dx'dx^ + H{p, z) dVL\_^, (4.1) 

and the coordinates = t, -(/', which includes spinning configurations as well (since the 
metric component g^t can be nonzero). 
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In what follows we show that, for d > 5, a rotating black hole with a spherical topology 
of the horizon can be written within this ansatz (then it corresponds to the Myers-Perry 
solution with a single angular momentum). Another interesting case is represented by a 
black hole with a 5^ x S'^~^ topology of the horizon rotating with respect to the azimuthal 
angle ip (thus with a rotating S"^). 

To make correspondence with (^5), it is convenient to choose the following parametriza- 



tion of (4.1 



ds^ = -Up, z)dt^ + Mp, z) {dp" + dz^) + hip, z){d^ + W{p, z)dtf + hip, z)dnl_^. (4.2) 

The resulting equations for /j, W still have a dependence on only two variables, their 
structure being quite similar to that found in the static case (and thus we shall not write 
them here). 

The rod structure as defined in Section 2 still holds for rotating solutions {e.g. /o(0, z) = 
still defines the position of a horizon). Moreover, for a generic configuration, the metric 
functions fi admit the same expansion at p = as in the static case. As p — )• 0, the new 
metric function W{p, z) has the following form: W{p, z) = wq{z) + p'^W2{z) + 0{p'^), with 
W(){z) = Qh = const, on a timelike rod, Qh being the event horizon velocity (note that the 
Killing vector d/dt + ^Ind/dip is null at the horizon). At infinity one imposes again the 
same asymptotic behaviour for fi, while W — t- 87rJ/Vrf_2(p^ + 2;^)^'^~^^''^, with J the angular 
momentum of the solutions. The mass of the solutions is read again from the asymptotic 
expression of /q. 

However, other kinds of black holes, such as the spinning Myers- Perry black holes with 
multiple angular momenta, and presumably also the spinning balanced black holes with 
5^ X S'^ horizon topology, do not fit within the conformal ansatz we use. For example, in 
the cases we studied, the conical sigularities of the spinning 5^ x 5^ black hole, could not 
be eliminated. This should not be taken as a sign of the absence of regular S'^ x S'^ black 
holes, but rather a consequence of the very restrictive form of the metric ansatz we employ. 

Another argument comes from the recently found new extremal near horizon geome- 
tries in Ref. [^^, that are similar to the near horizon geometries of the doubly equaly 
spinning extremal Myers-Perry black hole. These near horizon geometries were found as 
deformations of the doubly spinning extremal Myers-Perry black hole with equal spins and 
have angular cross terms that we do not consider in ( [4.1[ ). So, taking into account the sim- 
ilarities between these near horizon geometries, we expect that balanced rotating S"^ x 5^ 
black hole or even more exotic, less symmetric, rotating black hole solutions not to fall in 
the (|4.1|) category. 



4.1.1 The singly spinning Myers-Perry black hole 

One can show that a Myers- Perry black hole with a single nonvanishing angular momentum 



can be written with the line element (4T). This solution is usually written in the following 
form 1^] 

ds^ = -dt^ + {dt - a sin^ d^f + -dr^ + J^dO^ (4.3) 

+ (r^ + a^) sin^ 6 dtP^ + cos^ 9 dill_^ , 
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where 



S = + cos^ 9, A = + - , (4.4) 



such that at infinity the hne element (^]^) is approached. Employing the same change of 
coordinates ( p. 131 ) as in the static case, where now 

G{r) = 2j^dr, (4.5) 



we can get the metric to be in the form (4.1). Similar to the Schwarzschild-Tangherlini 
case, the above integral has a particularly simple form for d = 5 

G{r) = 21og 2 (^r + + - ^) , (4.6) 

and for d = 7, where 



G(r) = log + 2 [r'^ + Var^^ + r^^jJ • (4.7) 
Therefore the explicit form of the metric functions v, Gij and H in ( [4.1] ) can be derived 



from (|4.3| ). However, their expression is very complicated and we shall not present it here^*^. 

One may wonder whether more general Myers-Perry solutions can also be written 
within the ansatz (|4.1| ). For example, when n ^ N angular momenta are equal, the Myers- 
Perry black hole exhibits a symmetry enhancement to U{n) xU{l)x Rt symmetry. However, 
we have verified that such black holes require a more general metric form than ([4.2]). 



4.1.2 A rotating S"^: d = 6 black holes with S'^ x S'^ topology of the horizon 

The Myers-Perry black hole has a S"^^"^ horizon topology. However, in principle, solutions 
with a more complicated horizon topology can also be written within the ansatz ( [4.1D . 

For example, by employing the same methods as in the static case, we have constructed 
d = 6 rotating black holes with a S'^ x S"^ topology of the event horizon (note that only one 
of these two spheres has a round shape). They are found by starting from a static solution 
and increasing the value of the angular velocity iln of the event horizon, which enters the 
boundary conditions at /) = 0. Then, by varying also the second parameter a/b, associated 
with the position of the rods, the full set of 'generalized black rings' with a rotating can 
be explored in principle. As the second V'— rod extends to infinity {a/b — )■ 0), the radius of 
the horizon of the round S'^-sphere increases and asymptotically it becomes a two-plane. 
Here one expects to recover, after a suitable rescaling, the four dimensional Kerr black hole 
uplifted to six dimensions. 

Although we did not yet investigate the full parameter space, all solutions we have 
constructed so far possess a conical singularity on the finite ip-rod. This is not unexpected, 
since these configurations are natural higher dimensional counterparts of the d = 5 black 



case 

20 



Note that also the expression for p = 5 sin 2^ ^ '^^^a-li ^ j ^ ^ ~ \ cos 26 ^1 + ^rj;j ^re simple in this 
Their explicit form for d = 5 is given e.g. in 
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ring with a rotating 5^ found in [36|. Different from the balanced black ring in (which 
has g^pt = 0, (7,^4 7^ 0), the solution in ||3^ has a conical singularity inside the ring for any- 
allowed value of the angular momentum^^. Balanced black rings with a rotating 5^ exist 
only if they rotate along as well |35]. Therefore, we expect a similar result also for d > 5 
solutions. However, this class of configurations would have a rotating S'^~^ and then would 
not be described by the ansatz ([4.11). 

A description of the spinning d = 6 black holes with S'^ x S"^ topology of the horizon, 
together with generalizations for d = 7 and multi-black hole objects will be presented 
elsewhere. 

4.2 A rotating S"^""^ 

4.2.1 A general metric ansatz 

Heuristically, to provide a centrifugal force which may balance a system with a nonspherical 
horizon topology, one needs to rotate the S'^~^ sphere^^^. In principle, the simplest case 
corresponds to black objects with a single angular momentum with respect to a direction 
on S'^~'^. For d > 5, a possible generalization of the static ansatz ( p.8D to this case has a 
dependence on p, z and an angular variable 6 on S'^~'^, with 

ds2 ^ g^.(^p^ ^^i^^j + Q^.^p^ ^^i^^j + fj^p^ 2^ ^) dnl_Q, (4.8) 

where = p,z,6 and x' = t,il^,ip (note that the metric on 5"^~^ is written in terms of a 
warped product of S"^ and S'^~^, with dQ,'^_^ = dS^ + sin^ Odip"^ + cos'^ 9dQ,'^_Q in the absence 
of rotation). 

An explicit parametrization of the above line element which makes contact with the 
static ansatz (|2.8D, proposing also a choice of the gauge in the z^-sector, is 



ds^ = -Mp, z, e)dt^ + flip, z, 6) {dp^ + Slip, z, e)dz^) + f2{p, z, e)dil? (4.9) 
+ h{p, z, 9) (dO^ + S2ip, z, 9) sin^ e{d^ + W{p, z, e)dtf + S^^ip, z, 6) cos^ edal_ 



with Si = 1 and = in the static case (note that Ss = for d = 6; also we did not 
consider solutions with rotation on S"^ , i.e. g^t = 0). Therefore finding rotating solutions 
with a nonspherical topology of the horizon reduces to solving a set of coupled partial 
differential equations for fi,Si,W, with suitable boundary conditions. However, this is a 
very difficult numerical problem, since the equations depend on three variables. 

The only configurations we have attempted to construct within the above ansatz cor- 
respond to d = 6 black holes with S"^ x 5"^ topology of the horizon (this time both spheres 
deviate from sphericity). Again, one starts with static solutions and increases the value of 
the angular velocity ^h- On general grounds, balanced solutions are expected to exist for 
a critical value of ^h- 

However, although in principle all methods and the software used in the static case can 
also be applied for this 3D problem, so far we could not make any progress in this direction. 



^^The thermodynamical properties of this solution are examined in 



^For example, in d = 5 dimensions, the balanced ring is rotating with respect to the direction. 
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The main problem is that we could not assure the convergence of the numerical process in 
the presence of rotation, even for small values of fin- This problem may be related to the 
issue of the gauge choice in this case^^ . For the version (|4.9D of the generic ansatz ( [4. 81 ) , we 
have fixed the metric gauge by setting to zero the extradiagonal components gpz , gpo , Qze ■ 
This seems to be a natural generalization of the "conformal gauge" employed in ( |2.8D {i.e. 
gpp = 9zz and gpz = 0), which is the most convenient choice in a numerical approach. 
However, different from the static case discussed above, we could not prove the consistency 



of the proposed metric ansatz (O). In other words, for rotating solutions there is no 
obvious way to prove that the constrained equations are solved automatically via Bianchi 
identities (plus suitable boundary conditions), as for the case of solutions with dependence 
on p,z only (see the discussion in Section 2.1). This problem survives for other metric 
gauge choices we have considered. 

One should also remark that d > 5 balanced black rings with a single angular mo- 
mentum can be constructed in principle by using the ansatz ( |4.9| ). However, all technical 
issues pointed out above apply for ring solutions as well. Moreover, higher dimensional 
black rings can also have a rotating 5*^. Again, the exact solution is known only in five 
dimensions |^^. For d > 5, we expect such solutions to have non trivial cross term metric 
components on the S'^^^ part. With the (S"^~^-spins turned on, the solutions will presum- 
ably have R x [/(l)""^-*^ symmetry and be cohomogeneity d — (n + 2) (with n = [(d — 2)/2]). 
And, even when n angular momenta are set to be equal (and non-zero), the enhancement 
to a U{n) X U{1) rotational symmetry would not be enough to lead to a simple metric 
ansatz as it was the case for the Myers-Perry black hole with equal angular momenta (see 
e.g. 0). 

4.2.2 d = 7 rotating solutions with S"^ x horizon topology 

The only possibility we have found so far to construct nonperturbative spinning black 
objects with rotation on S'^~'^ and possessing a nontrivial topology of the horizon, corre- 
sponds to the case d — 4: = 2k + 1 (with k = 1,2,...). There the problem can be greatly 
simplified, when the a priori independent {d — 3) /2 angular momenta on 5°^"^ are chosen 



to have equal magnitude, since this factorizes the angular dependence ||371]. The problem 
then reduces to studying the solutions of a set of five partial differential equations with 
dependence only on the variables p, z. 

In what follows, we present some partial results for the simplest case d = 7. The metric 
ansatz in this case is a straightforward reduction of (|4.g|), with 



ds^ = -Up, z)df + flip, z){dp' + dz^) + h{p, z)d^^ + h{p, z)de^ (4.10) 
+ U{P, z) ( sin^ 0{dipi - W{p, z)dtf + cos^ e{dip2 - W{p, z)dtf) 
- ifiiP, z) - h{p,z))svo? 6 cos^ 6 {difi - dip2f, 

where < < vr/2, < ^i,ip2 < 27r. The static ansatz ( p.^ ) is recovered for W = and 
/4 = /s- 



^^Note also that the complexity of the equations increases tremendously as compared to the static limit 

(El, (11). 
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A suitable combination of the Einstein equations, G\ = 0, Gp + Gl = 0, Gq = 0, 
= 0, = and = 0, yields the following set of equations for the metric 

functions: 

VVo - :^{Vfof + ^(Vfo) ■ (V/2) + ^(V/o) • (V/3) + Tr^(V/o) • (V/4) - hiVWf = 0, 

VVi - |(V/i)2 - -^(V/3)2 - ^{VWf - /^(V/o) • (V/3) - /^(V/3) • (V/4) 
/l . /l , /l 



(V/2) • (V/3) - TrjV(V/o) • (V/2) - Tr4V(V/o) • (V/4) 



/2/3' V 2/0/2^ ' 2/0/ 

/i .V7.^ .V7.^ , 2/! 



2 



2/2/4^^^'^ • (V/4) + -^(4/3 - /4) = 0, 

^'•^2 - ^(V/2)' + ^(V/o) • (V/2) + -^(V/2) • (V/3) + ^(V/2) • (V/4) = 0, (4.11) 
VV3 + ^{^k) ■ (V/3) + T^(V/2) • (V/3) + -^(V/3) • (V/4) + 4/i(^ - 2) = 0, 

^/O ^/2 ^J4 J3 

V'/4 - ;^(V/4)=^ + ^{VWf + -V(V/o) • (V/4) + ;rV(V/2) • (V/4) 



2/4^ fo' ' 2/0^ ^ 2/2 

|(V/3)-(V/4)-^ 
73 J3 



v'M^ - 7T^(v/o) • (viy) + -^(v/2) • (viy) + ^(v/3) • {vw) + -^(v/4) • {vw) = 0. 

^/O ^J2 /3 ^J4 

All other Einstein equations except for G^ = and Gp — G| = are linear combinations 
of those used to derive the above equations or are identically zero. Moreover, a similar 
reasoning to that presented for static solutions assures that the constraints G^ = and 
Gp — G^ = vanish identically via Bianchi identities plus suitable boundary conditions. 

The only rotating solutions we have attempted to construct within this ansatz have 
an S"^ X 5'^ topology of the horizon (i.e. d = 7 generalized black rings) and a rod structure 
similar to that discussed in Section 2 {e.g. the horizon is located at p = and —a < z < a, 
where /o = 0). Moreover, the expansion at p = of the metric functions fi (i = 0, . . . , 3) is 
similar to that presented in the static case. For p — )• and z < b, the new functions /4, W 
have the following expansion 

hip, Z) = Uoiz) + p'/42(^) + //44(^) + . . . , W{p, Z) = Wq{z) + p'w2{z) + . . . ,(4.12) 

where Wq{z) = const, for —a < z < a (i.e. on the event horizon). The expansion of 
f4,{p,z) is different for the il— rod (i.e. p = and a z-interval [6,00]), where fi{p,z) = 
p^/io(z) + p'^fiiiz) + .... The obvious boundary conditions for large p,z are that fi 
approach the Minkowski background functions ( 2.26D , while W = 0. 

The Killing vector x = d/dt + Qid/dipi + ^2d/d'-P2 is orthogonal to and null on the 
horizon. For the solutions within the ansatz ( |2.8| ), the event horizon angular velocities are 
equal, ^1 = ^2 = W{Q,z)\_a<z<a = ^H- 

As in the case of Myers-Perry black holes, these rotating black holes have an ergosurface 
inside of which observers cannot remain stationary, and will move in the direction of the 
rotation. 
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Figure 9. The metric function W (left) and gu — —fo + fiW^ (right) are shown close to the 
horizon for a d = 7 rotating solution with 5^ x event horizon topology. The input parameters 
here are a — 1, b = 4 and fin = 0.05. The contour line in the right panel indicates the ergo-region. 



The ergosurface is located at gu = 0, i.e. 



/o + fiW^ = 0, 



(4.13) 



and does not intersect the horizon. 

The area Ah and the Hawking temperature of the black hole can be expressed as 



Ah = ^Ai^ dzh./hhh 



Th 



lim 



fo{p,z) 



p=0 



2tt p^o W p^fi{p,z) 



(4.14) 



-a<z<a 



The ADM mass M and the angular momenta Ji = J2 = J of the solutions can be read 
from the asymptotic expression of fo,W: 



/o ~ 1 



levrM 



+ ..., W 



SvrJ 



These spinning black objects satisfy the Smarr formula 

4 1 5 

-M = -ThAh + -^hJ . 

5 4 2 



(4.15) 



(4.16) 



The rotating solutions are found by starting with the static configurations with S"^ x 5^ 
topology of the horizon discussed in Section [j.1.2 and increasing the horizon velocity ^}h- 
We have found that the absolute value of the conical excess 5 decreases with Qh, which 
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suggests the existence of a critical value of the horizon angular velocity such that 5 = (this 
is the case for the general d = 5 Emparan-Reall rotating black ring |18|, see Appendix B). 
Unfortunately, the accuracy is lost and the numerical process diverges before approaching a 
balanced configuration, without being possible to identify a clear origin of this behaviour. 
However, a more sophisticated numerical approach may be able to find balanced black 
holes with x topology of the horizon. 

As an example, we show in Figure 9 the metric function W{p,z) and gtt{p,z) for the 
rotating solution with parameters a = 1, 6 = 4 and fin = 0.05 (the shape of the other 
metric functions is similar to that found in the static case). In the right panel, one notices 
the existence of a region in the {p, z) plane with gu < and of an ergosurface where 

/O = f4W\ 



5. Conclusions 

The main purpose of this work was to present a general framework for the nonperturbative 
construction of a class of d > 5 static black objects with a nonspherical topology of the 
horizon. The solutions are found by solving numerically a set of four partial differential 
equations with suitable boundary conditions. Such an approach may be viewed as comple- 
mentary to the approximate construction of such black objects developed recently in Q, 
|1C], since it may work well if the length scales involved are not widely separated. Also, 



this made possible to consider some black object topologies that are not captured within 
the blackring/fold approach [e.g. S"^ x S'^ in d = 6). 

As a concrete application of the proposed formalism, we have presented numerical 
evidence for the existence of several d > 5 black objects with a nonspherical topology of 
the horizon. These solutions represent generalizations of the d = 5 static black rings, dirings 
and Saturns, with similar basic properties. Without entering into details, we mention that 
the double analytic continuation •i/^ — )• iT, t ^ ir in the line element (|2.8| ) leads to the 
more exotic interpretation of the solutions in this work as bubble-black hole sequences in a 
Kaluza-Klein theory. For example, the black hole with a S"^ x S'^~^ topology of the horizon 
becomes a pair of black objects (with one accelerated horizon) sitting on a bubble. 

Not completely unexpected, our static solutions always possess conical singularities. 
The only way to achieve balance seems to be to rotate the solutions, no other mechanism 



being known at this moment. For example, the arguments in [12] put forward for general- 
ized black rings apply directly to all static solutions in this work and one can show that the 
conical singularities plague also the Einstein-Maxwell-dilaton generalizations of the solu- 



tions in this paper^^. Moreover, the results in [39| show that the Gauss-Bonnet corrections 
to Einstein gravity cannot eliminate the conical singularity of a d = 5 static black ring, and 
we expect a similar result to hold also for the higher dimensional configurations discussed 
in this work. 



^^Static balanced black objects with a nonspherical topology of the horizon may exist, however, if the 
gauge fields are not vanishing at infinity, i.e. for an asymptotic Melvin structure of spacetime Such 
solutions are known in closed form in d = 5 dimensions, see e.g. |Bq|. 
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As argued in Section 4, the construction of the rotating balanced version of the solutions 
we have considered is a much more difficult task. However, based on the experience with d = 
5 exact solutions, it is likely that some of the qualitative features of the static configurations 
will hold also in the spinning case. 

A further generalization of the solutions may be along the lines of Ref. |40|, where 
(apparent) horizons of topology x 5""+^, n, m > 1 were considered. 

Thus, we expect that the new configurations discussed in this work represent just 'the 
tip of the iceberg' and a variety of new d > 5 black objects with nonspherical topology of the 
horizon are likely to be discovered within a nonperturbative approach. In any such attempt, 
the rod structure of the solutions (or a suitable generalization of it) would represent an 
important ingredient, as a tool to fix the topology of the horizon. For example, it would be 
interesting to adapt the numerical methods in this work for the domain structure approach 
introduced recently in [17|. 

In our opinion, any progress in this direction would require the development of a 
consistent numerical scheme capable to solve as a boundary value problem the Einstein 
equations with a dependence on at least three coordinates. 
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A. Five-dimensional seeds 

For completeness and comparison with the higher dimensional counterparts, we present in 
what follows the expression (in a form suitable for numerical calculations) and some basic 
properties of the five dimensional seed solutions^^. Also, the functions fi below are used 
as background functions for the corresponding higher dimensional solutions. 

Moreover, based on the results below, one can easily construct e.g. the d = 5 counter- 
parts of the Figures 6 and 8, which clearly show that the d = 5 pattern repeats in higher 
dimensions. 

A.l The static black ring 

The metric functions fi of the static black ring are given by [P,|18] 



f _ P2+2g2 f _ (Pl+2gi+P2)(Pl+P2+6+6+c(-gl+6+2g3-Pl+P2+2P3) /a ^ 

^0 - Pi+25l ) Jl - 8(l+c)(Pl+Ci)(P2+52)(P3+?3) ' ^ 

/2 = §(^^3 + 26), /3 = ^'3, 



■^^To the best of our knowledge, this study is missing in the literature and thus may be useful for future 
studies. 
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where = z — Zi, 



Pi = ^/p^ + {z-z.i)^-{z-Zi), (A.2) 

and 

zi = -a, Z2 = a, Z3 = b, (A. 3) 

a and b being two positive constants, with a < b and c = a/b. 
The leading order expansion as p — >• of these functions is: 

„ z + a a-z p'^{z-a) 
JO — , Ji — 77 — ; — /: 



z-a 2{b-z){a + z)' 2{b-z){a + z)' 

/a = 2(6 — z), for — 00 < 2; < —a, 

p2 2(a^ — z^) 

4(a2-z2)' " (a + 6)(a2-z2)' " (6 - z) ' ^^'^^ 

/3 = 2(6 — z), for — a < z < a, 

, , 2a (6-a)(z + a) /^^(z + a) 

/O — -L ; — , Jl = TTTT"; — T7 \7I 7' /- ~ 



z + a 2(6 + a)(z-a)(6-z)' 2(z-a)(6-z)' 

/s = 2(6 -z), for a < z < 6, 
„ 1 2a „ 1 „ 2(z-6)(z + a) /)2 . 

z + a 2(z — 0) (z — a) 2(z — 6) 

The mass, event horizon area and temperature of the d = 5 static black ring are: 



These black rings have a conical deficit for the finite ip-rod, with 



A.2 The static Saturn 

The metric functions are given in this case by 

r _ P1P3 . {P^ + PiP2)Hp^ + P2P3)Hp^ + + P3A)A , . 

P2PA ' (p2 + p2)(^2 + p2)(^2 + p2)(^2 + p2)(p2 + ^,^3)2(^2 + p^p^^ ' ^ ' ^ 

f - ^^^2 f _ p 



where Pi is given by ( [A .21 ) and 

zi = —a, Z2 = a, Z3 = 6, Z4 = c, (A. 8) 

a, 6 and c being three positive constant, with a < 6 < c. The notation here is somehow 



arbitrary and has been chosen to make contact with the results in |12 
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This describes a multi-black hole solution, with a black ring with horizon topology 
S"^ X around an black hole. 

The leading order expansion as p ^ of these functions is: 

{z + a){z-b) z-a _ - a) 

Jo = 7 w /I — T77I V? — : — \ ' r- 



{a-z)ic-z)' ' 2{b-z){a + z)' ' 2(6-z)(a + z)' 

fs = 2(c — z), for — oo < z < —a, 

p'jh-z) 2a''{a + c){h-z) _ 2(a^ - z'') 

4(c-^)(a2-z2)' •^i (a + 6)2(c-z)(a2-z2)' •'2 _ ^) ' 

/3 = 2(c — z), for — a < 2; < a, 

^ {a-z){z-b) ^ {a-bf{a + c){a + z) 

(c-z)(z + a)' ■'^ 2{a + b)^a-c){a- z){b~ z)' ^ ^ 

/2 = ./'^^t"^ /3 = 2(c-^), fora<z<6, 
2(a — — b) 

{a-z)p^ {a + c){a-z) 

JO — -TTi 77 77 : 7' h 



A{b- z){c- z){z + ay ' 2{a - c){b - z){a + zY 

/2 = ^^^^^^,/, = 2(c-.),tor6<.<c, 

{z-a) 

. {a-z){c-z) „ 1 , 2{b-z)iz + a) , r ^ 

/o = ^ r(, /i = 7, /2 = ' h = 0/ N > for c < z < 00. 

(2; + a)(2; — 0) 2(z — c) (a ~ -2^) 2{z — c) 

The event horizon area and temperature of the black hole with horizon topology S'^ x 
are: 



^g)^^^^,V2aVa±l (A.10) 
The same quantities for the black hole with horizon topology are: 

^ V c-a ' \/2^/27r V (a + c)(c- 6) ' 

These solutions have a conical deficit for the finite '0-rod {i.e. p = 0, —a<z<a), with 



, b + a c — a\ . . X 



which prevents the configuration from collapsing. 

The ADM mass of this system, as measured at infinity is: 



M = —{2a-b + c). (A.13) 



A. 3 The static di-ring 

The metric functions are given in this case by 



/0=^, /2 = ^, /3 = P5, (A.14) 



P2P4' •'^ ~ P1P3P5 

P5(p'+PlP2)'(p'+PlP3f(p'+PlP4)"(p"+i^'iP4)"(p"+P2P5)(p'+P4/^) 
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where Pi is still given by ([AT2|) and 



zi = —d, Z2 = — c, 2:3 = —a, Z4 = a, Z5 = 6, (A. 15) 

a, 6, c and d being three positive constants, with d > c > a and a < b. This describes a 
configuration consisting of two concentric black rings with horizon topology S'^ x S^. 
The leading order expansion as p — )• of the metric functions is: 

_{z + a){d + z) _ {a-z){c + z) 

~ (z-a)(c + z)' ~ 2(a + z)(z-6)(d + z)' 

^2 = ' w A ' /3 = 2 6 - z , for - 00 < z < -d, 

2{a + z){z — b){d + z) 

_ p^{a + z) _ {c-df{a + df {a + z) 

~ 4(a -z){c + z){d + zy 2{a - df{h + d) {a - z){c + + z) ' 

2{a - z){c + z){d + z) J/ ^ 

/2 = 77 w — \ — ^ , h = 2 6 - 2; , for -d<z< -c, 

(5 - z)(a + z) 

_ {a + z){c + z) (a + d)2 (a + 2;)(rf + z) 

~ (z - a)(z + d) ' ~ 2(a - d)2(6 + d) (a - z)(c + z) ' 

/2 = TTT — -^^7 u\( I ^ ' -^3 = 2 6 - z , for -c<z< -a, 

2(a + — o)(c + z) 

{c + z)p^ {b + c){a + d)^ {a-z){c + z) 

n = -K} , , A 7 ^ — T7 rTTTT — T' (A.16) 



4(d + z)(a2 -z2)' •'^ 2(a + c)2(6 + d) (a + z)(z-6)(d + z)' 
2(a2-z2)(d + z) „ . , 

(6-z)(c + z) ' ^3 = 2(6-.), for - a < z < a, 



_ {c + z){z-a) _ {b + c){a + z){d + z) 
{z + a)iz + d) ' ~ 2(a + + - a)(c + z) ' 

•^2 = w ,w ' V -^3 = 2 6 -z), for a < z < 6, 

2(a — — o)(c + z) 

^ {c + z){z-a) ^ 1 

{z + a){z + d)' 2(z-6)' 

, 2{a + z){z -b){d + z) „ p2 

^2 = ( u \_ \ ^ ^3 = for 6 < z < cx). 

(z — a)(c + z) 2(z — 0) 

The event horizon area and temperature of the "left" black ring (with the horizon located 
at p = 0, —d < z < — c) are: 

^(5) ^ ^^2 V2(c-d)2(a + d) 1 ^(5) ^ 1 d-a^b + d 
^ (d-a) ^/bTd' ^ 2V27rd + a d-c ' 

The same quantities for the second black ring are: 



(a + c) Y(a + 6)(6 + d)' 4\/27r a(a + d) V 6 + c 
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These solutions have a conical deficit for the third -i/^-rod with a < z < b, with 



which prevents the configuration from collapsing. 

The ADM mass of this system, as measured at infinity is: 

Svr 

M = —{2a-c + d). (A.20) 

B. 'Isotropic' coordinates and new diagrams 
B.l A coordinate transformation 

The coordinates (p, z) defined in Section 2 have the advantage to make contact with d = 5 
Weyl coordinates and to make possible to visualise some basic properties of the solutions 
in terms of rod diagrams. However, although the domain of integration has a rectangular 
shape, the range of both p and z is unlimited and, within our numerical scheme, it is rather 
difficult to construct suitable meshes, especially in the z— direction. 

In practice, we have found another coordinate system which has proven useful in 
the construction of some black holes with a nonspherical topology of the horizon. The 
transformation between (p, z) and the new coordinates (r, 0) goes as follows. Starting with 
one of the diagrams in Figure 1, let us choose a finite rod there as 'central rod' {i.e. that 
extends from —u < z <u). Then we introduce the coordinate transformation 

p{r,e) = ^sm2e, z{r,e) = ^-^cos2e, (B.l) 

with Tq = n and tq < r < oo, < < ^, such that the 'central rod' is located at r = ro 
and all 0-interval. 

Then all static solutions in this work (including the d = 5 configurations^^) can be 
studied with a metric ansatz akin to (|2.8D , with 

ds^ = -fo{r,e)dt'' + fi{r,e){dr^ + r^de^) + f2{r,e)dij^ + Mr,e)dnl_^. (B.2) 

For ro = (and only two semi- infinite rods), one recovers the flat spacetime metric with 
/o = /i = 1, /2 = cos^ 6 and /s = cos^ 9. The expression of fi for a more complicated 
rod structure can easily be derived once we know the solutions in (p, z) coordinates (note 
that /i in (p.2|) does not coincide with fi in (2.8), since the Jacobian of the transformation 



(B.l) enters there also). 

Also, it may be interesting to remark that r and can be viewed as 'generalized 
isotropic coordinates'^^. This is justified by the observation that for the simplest case of 



^This approach can be extended to d = 4 axisymmetric solutions. However, in this case the transforma- 



tion between the Weyl coordinates p, z and the spherical coordinates r, 6 is not given by (B.l). 

^^It may be interesting to notice that a version of the isotropic coordinates has been used in most of the 
previous numerical studies on d = 4 asymptotically flat axisymmetric solutions, see e.g. [B0|. 
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a single black hole with spherical topology of the horizon, the Schwarzschild-Tangherlini 
solution in isotropic coordinates is recovered (the results for d = 5 below can easily be 
generalized to higher dimensions). 

One of the advantages in the numerics of these coordinates is that the range of 6 is 
finite. Moreover, the coordinate singularities in the metric functions are easier to handle in 
this case. For example, we have constructed in this way the full set of d = 6, 7 black objects 
with a 5^ X S'^~'^ topology of the horizon, with a better accuracy than that obtained for 
the {p, z) coordinate system. 

The new coordinate system leads also to a new type of diagrams, which is the coun- 
terpart in (r, ^) coordinates of the rod-diagrams in Figure 1. As one can see from (|B.lJ ), 
p = corresponds to r = tq or = 0, 7r/2. This suggests to show the domain of integration 
together with the boundary conditions satisfied by the metric functions /o,/2 and /s. In 
our conventions, a wavy line indicates a horizon /o = 0, a thick line means /2 = (i.e. a 
V'— rod) and a double thin line stands for an 0— rod, /a = 0, see Figure 10 (the general- 
ization of the d = 5 diagrams there to higher dimensions is straightforward). The horizon 
topology can also easily be read from that figure: a spherical horizon continues with rods 
of different directions, while for a black ring, the horizon continues with -;/'— rods only. 

These features are clearly illustrated by a number of d = 5 exact solutions which we 
shall present in what follows. 



B.2 d = 5 static solutions in 'isotropic' coordinates 

We shall start with the simplest example, corresponding to a Schwarzschild-Tangherlini 
black hole. The metric functions in this case read 



,r 



/o(0 = (^)^ /iW = ^i + 3j , h{r,e) = [i+'f,] cos'e, (b.3) 



2\ 2 



„2\ 2 



f3ir,0) 



.2\ 2 



1 + 



sm 



with an event horizon at r = rg, the corresponding diagram being shown in Figure 10a. 

The Emparan-Reall static black ring has also a relatively simple expression in these 
coordinates. 



fair) = ( 



r — r. 



+ r. 



^2 
2> ' 



f2{r,e) 



r"^ + rQ\'^ sin^ 6 cos^ ' 



h{r,e) 



/3(r,^) = - 2i?3+r,^ 1+ - 



— 1 cos 20 



r / V ro 



(B.4) 



r- / 1 + ^ 
-Ks \ 1 + ^ 



where 



= - 



b 



l+(- 

r 



(1 + (-)')(! + (-)') + 2(^)'(^) - 2(^)2 cos 20 



„2 \ 4 



2 \ 2^ 



2cos20f^V) il+l^y -2cos2e(^) 
r ) ) \ \rn) \rn) 



1/2 
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Figure 10. The domain of integration for the coordinate system (B.2) is shown for a Schwarzschild 
black hole, a static black ring, a black Saturn and a diring in d = 5 dimensions. 



(Note that the metric function gu has no angular dependence.) The new type of rod- 
diagram for a static black ring is shown in Figure 10b. The horizon is again located at 
r = ro, the finite ip-iod with an angular excess being at 6 = 0, vq < r < rt,. 

The black Saturn can also be written in (r, 0)-coordinates, with the following expression 
of the metric functions 

(4^)'!^' /2(r-,g)= ^''V°^' |^cos^g, /3(r,g)= ^"'V°^' sin'e, (B.5) 

/ / R2 + ^^^^ cos^ g \ ' Ri + ^^^^ cos^ 6 {p^{r, 9) + R^R^Y 

V ^r^J Ui + ^^cos^^J i?2 + ^^cos2e(p'(^>^) + ^?)(p'(^>^) + ^i)' 



/o(r,0) = 

h{r,e) = 
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where we have introduced the auxihary functions 

R,{r, 9) = { - + r^) - rl{r^ + 4) cos 26 



+ \l {r'^ + ri + 2r2r2 cos 26*) (rg + r-^rf + 2r2r^r2 cos 26) 



R2{r, 6) = ^ - r\rl + r^) - rj(r^ + rl) cos 26 

+ {r'^ + + 2r2r2 cos 26*) (rg + r'^r^ + 2r2r^r2 cos 26) 

and p{r,6) as given by ( p.l| ). The new rod-diagram for a static black Saturn is shown 
in Figure 10c. One can notice the existence of two horizons, at r = rg and at 6 = vr/2, 
i"d ^ r < rs, respectively. 

Finally, we give also the expression of the metric functions for a static di-ring in (r, 6)- 
isotropic coordinates (the corresponding diagram is shown in Figure lOd) 



'R^ + ^^^S^ cos2 e\ {R^R, + p2(^^0))(i?^i?4 + p^(r,6))^ 



\ Rl + p^r,e) J {Rl + p^{r,e)){Rl + p^{r,6){RiR5 + pHr,6)){RsR5 + p^{r,6))' 
The auxiliary functions Ri have the following expression 

^) = ^ ( - ^'(^0 + ) - + 4)r^s cos 20 



+ a/ (r"^ + + 2r2r2 cos 26*) (rg + r^rf + 2r'^r^r'^ cos 26*) 



0) = ^ ( r\4 + r^) - (r^ + 4)rj cos 20 



1 



+ ^J (r^ + - 2r2r^ cos 29) (rg + r^r^ - 2r'^rjr'^ cos 20) j , 

ii5(r, 0) = (^r2(r^ + rj) - (r^ + 4)rj cos 20 

+ ^(r4 + rj- 2r2r2 cos 20) (rg + r'^rj - 2r'^rjr^ cos 20)^ , 

Different from the previous cases, the line r = tq corresponds in this case to a finite ip-rod 
with a conical excess. The first event horizon is located at = 0, tq < r < r^, and the 
second one at = 7r/2, vq < r < r^. 

These examples make clear that the functions fi have a manageable expression also in 
(r, 0)-coordinates. In particular, the singularities in the expression of /i at p = 0, z = ±it 
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are eliminated by the coordinate transformation ( [B.lD . Also, it is straightforward to per- 
form a systematic study of these solutions in these coordinates, similar to that considered 
in Appendix A. 

B.3 Rotating solutions: balanced black ring and Myers-Perry black hole 

The coordinate system introduced above turns out to be very useful also in the numerical 
construction of d = 5 rotating solutions. We illustrate that by exhibiting some results for 
the balanced Emparan-Reall black rings and Myers- Perry black holes with a single angular 
momentum. 

The spinning solutions can be constructed within a simple generalization of (|B.2|) , with 



ds^ = -Mr,e)dt^ + T7^(dr2 ^ ^2^^2) ^ /^(r,^)^^' + f3{r,9){dv + ^^^dtf, (B.7) 
Ji{r,0) r 

such that the horizon is located at a fixed value of r = tq. Expanding the Einstein equations 
in the vicinity of the horizon in powers of r — r^, one finds fi{r, 9) = fiQ{0) + fi2{0){r — rh)^ + 
0{r — r/i)^, w{r,6) = Wh + W2{0){r — rh)"^ + 0{r — rh)^, (where the functions fik{9),W2{9) 
are solutions of a complicated set of nonlinear second order ordinary differential equations 
and fooiO) = 0), which leads to an event horizon metric 

^2 j/)2 

= T^m + /2o(^)#' + /3o(e) V- (B.8) 
Jl0[9) 

The Hawking temperature, entropy and the event horizon velocity^^ of the solutions are 
given by 



Th = ^VfMKm, s = vrVo rdeJl^^fM^, = (b.9) 

2vr Je=o y fio{0) vq 

For any topology of the horizon, as r — )• oo, the Minkowski spacetime background is 
recovered, with /o = /i = 1, /2 = r^cos^^, /s = w = 0. The mass M and the 

angular momentum J of the solutions are read from the asymptotic expansion of the metric 
functions, /o = 1 - SttM /S-rrr^ + w = AJ /nr"^ + .... 

However, the expression of the metric functions fi,w are quite complicated for any 
topology of the horizon. For example, a straightforward but cumbersome computation 



based on the Weyl-coordinate expressions in |18] leads to the following metric functions of 



a balanced black ring in the (r, 0)-coordinates proposed above: 

JO - , ^2^2 TTT^ — ' - - '^o) r /2 - -7TI2I27I2 ^^^^ (^-l") 



r, 



0) 



(r2 



2-r2)2 Q _ r^onir^ +ro?iro+r^b)\f^{+^ S2S3 



2r2r2(r2+r2)2 5i53' {r^ - rl? {rl - rl) Q 



^*Note that the KiUing vector d/dt + Qnd/dLp is null at the horizon. Also, the Einstein equation = 
implies that the Hawking temperature is constant. 
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where, in order to simplify the expresssions, we have defined 

51 = (r^ + 4rVg + r^)(r^ + r^) + ^r'^ry^ cos 26 - ^r'^rlrlR:^, 

52 = -{r^ + 4)r! + r\r^ + r^f - 2r\y cos 20 + 2ryR3, 
^3 = (r^ - 4rVg + r^)(4 + 4) + 4rV^r^ cos 20 - 4rVgr^i?3, 

^4 = {rt + - rl{r^ + r^) cos 29 + 2r^V2i?3, (B.ll) 
^5 = (r^ + ^r^rl + r^)(r^ + r^) + 4r2r^rfe2 3^ + 4rVgrb2i?3, 
>56 = iri + r^)r2 - rl{r'' + r^) cos 20 - 2r2r2i?3, 



and 



^7 = {r" + r;^)(r^ + r^) - ^r'rlrt cos 20 + ^ry^Rz, 



Ui = (r^ + rl){rl - rtf + 2r2r2(3r^ + rt){rl + Sr^^) + IQr^rlrUrl + r^^) cos 20, 



,2 c. 16r-V^ (r^ + rg)'(rg + rg)' (rg + r^) 
(r2 -r2)2(r2 -r2] 



n—C'^Q OV ^ ' 0/' V ^ ' fc/' v I ' c2 m 10\ 

^/ — D]^04 — ^ 2\2/ 2 2\2 2' l^-J^^/ 



i?3 = _J_ J (r4 + - 2r2r2 cos 20) (rg + r^r^ - 2r'^ry cos 20) . 
2r 

The physical quantities are complicated functions of the input parameters ro,rft 
M = 37r^^^^^, J = V2tt — , Qh = — , ,(B.13) 



(^0 - n? n{r^ - nf ^(^2 



2^2 



The existence of this exact solution allows us to test the scheme developed in this 
work by recovering numerically the balanced black ring starting with the static solution. 
Then the Einstein equations are solved for the metric ansatz ( |B.7| ) and the rod structure in 
Figure 10b, looking for balanced solutions. Again, in practice we use a set of background 
functions which takes automatically into account the sets of conditions on the boundaries 
and determines the topology of the horizon. Therefore one defines fi = Fiff'\ where 
ff'^ are the functions of the static black ring as given by (|B.4|). In our approach, the 



position of the horizon tq and the radius of the ring r;, are kept fixed and one varies the 
event horizon velocity VLh- All other relevant quantities are evaluated from the numerical 
output. When increasing the boundary parameter Wh-, the absolute value of the angular 
deficit excess decreases, such that 5 becomes zero for a critical value of the event horizon 
velocity (afterwards the ring becomes over-rotating with 5 > 0). By varying the value of 
(or the position of the horizon), the full spectrum of Emparan-Reall balanced black rings 
can be recovered numerically. 

This approach turns out to be consistent and usually provides very good accuracy re- 
sults2^ . A crucial ingridient of our approach is that all numerical singularities are absorbed 

^^Note, however, that the black rings with a large radius or those close to the naked singularity point are 
difficult to obtain with enough accuracy. 
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Figure 11 (continued). 



already by the background functions of the static solution, such that the rotation leads 
to smooth functions Fi,w. The description of the numerical method presented in Section 
2 is also valid in this case. In particular, we have used a compactified radial coordinate 
X = 1 — ro/r and a nonequidistant grid in 6. 

In Figure 11, we plot the metric functions of an exact solution with tq = 1, = 1.93, 
as well as the difference between the exact solution and the numerical result. One can see 
that the differences are on the order of 10~^ everywhere. As shown in Table 1, the global 



quantities computed numerically according to (B^) are also in excellent agreement with 
the theory values ( B.13D . 



For completeness, we give here the expression of the metric functions which enter the 



Myers- Perry solution with a single angular momentum within the coordinate system (B.7 



.2 



/o = (l-(^m /i = ((l + (^)Y + ^cos2 0r\ h = r\os'0il + C-^ff, f3 = r'sin'eU, 



^ = r(z^ + -3-)-FT^ and [/ = (1 + (-) ) + + . " . . — , (B.14) 



[/' "-r'' r2 r2 (1 + (10)2)2 + a^cog2 0' 



the relevant quantities being given by 



M = ^a(a2 + 4r2), J = ^^a{a^ + Arl) , J]^ = (B.15) 



We have verified that within the same numerical scheme as above (with the background 
functions ff'^ given by the expressions (|B.3| ) of the Schwarzschild-Tangherlini metric in 
isotropic coordinates) one recovers the set of rotating black holes with an topology of 
the horizon and a single angular momentum. The numerical accuracy is even better in this 
case, since the expression of the background functions contains no square roots. 

It would be interesting to recover within the same approach the balanced black Saturn 
and balanced black diring solutions starting with the corresponding static configurations. 

To conclude, we have proposed a numerical scheme which could reproduce physically 
interesting d = 5 spinning solutions starting with the corresponding static configurations. 
This opens the possibility to study generalizations of the Emparan-Reall balanced black 
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109.545 
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1.93186 


0.204124 


16.0000 


16.0000 
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58.7878 
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446.645 


2.18890 


0.207020 


13.3332 


13.3333 


45.0836 


45.0843 


332.911 


332.909 
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n on /1 1 0/( 


IZ.UUUl 


Iz.OUOU 
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39.1918 


2i6.bL4: 


zYo.olo 


2.80588 


0.193649 


10.6671 


10.6667 


34.4289 


34.4265 


210.552 


210.563 


3.14626 


0.182574 


9.99982 


10.0000 


32.8624 


32.8634 


176.553 


176.555 


3.45197 


0.172516 


9.59981 


9.60000 


32.4596 


32.4607 


154.723 


154.726 


3.99215 


0.155902 


9.14274 


9.14286 


32.9869 


32.9877 


127.614 


127.616 


4.46653 


0.143019 


8.88879 


8.88889 


34.1828 


34.1834 


110.970 


110.973 



Table 1. The values of the event horizon velocity VLr, mass parameter M, angular momentum J 
and of the event horizon area Ah are shown for rotating balanced black ring solutions with tq = 1 
and several values of r^,. For comparison, both the numerical and exact values are given here. 

rings and Myers-Perry black holes in various theories where closed form solutions are un- 
likely to exist {e.g. in Einstcin-Gauss-Bonnet or Einstein- Yang- Mills theory). We hope to 
report on that in future work. 
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